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We consider an alternative warm inflationary scenario in whicii n scalar fields coupled to a dis- 
sipative matter fluid cooperate to produce power-law inflation. The scalar flelds are driven by an 
exponential potential and the bulk dissipative pressure coefBcient is linear in the expansion rate. We 
find that the entropy of the fluid attains its asymptotic value in a characteristic time proportional 
to the square of the number of fields. This scenario remains nearly isothermal along the infiationary 
stage. The perturbations in energy density and entropy are studied in the long-wavelength regime 
C . and seen to grow roughly as the square of the scale factor. They are shown to be compatible with 

COBE measurements of the fiuctuations in temperature of the CMB. 

I. INTRODUCTION 

> , 

CS| I The very early Universe was supposedly populated by a host of scalar fields but soon only one became to dominate 
' the dynamics, the others settled in the minimum of its potential. Of particular interest from the point of view of 
cosmological inflation are scalar fields with exponential potentials since these are natural candidates to drive power- 
LiJ , law inflation in Friedmann-Lemaitre-Robertson- Walker (FLRW) universes, i.e., a{t) oc where a{t) is the scale 
factor and a ~ constant > 1 [Q] . 
, As is well known scalar fields possessing exponential potentials appear naturally in different theories of fundamental 
• physics as superstrings and higher dimensional theories [H , in TV = 2 supergravity ^ as well as in theories undergoing 
"■5 \ dimensional reduction to an effective four-dimensional theory Q|. However, in many cases they happen to be too 
II • steep and fail to do the job as a < 1. Nevertheless, it has been shown that if one considers n non-interacting scalar 
O ' fields with exponential potentials they can cooperate to achieve power-law inflation in spite of the fact that no single 
^ . field can achieve it by itself The rationale behind it is that while each field descends toward the minimum of its 
^ ' potential the cosmic expansion rate (to which all the fields cooperate) acts a friction force upon each of them. This 
. . , result has been extended to Bianchi I and VIIq cosmologies as well |6| . A further interesting feature is that the larger 
^ ' n, the closer the resulting spectrum of initial cosmic perturbations is to Harrison-Zeldovich's. The proponents of this 
k>( \ scenario termed it "assisted inflation" ; however we find more fitting to call it "synergistic inflation" , and we shall do 
■ so henceforth. 

I A more realistic problem arises when in addition to the scalar fields one considers a matter fluid. The price to be 
paid is a rather involved set of field equations even in the simplest case when the fluid interacts with the scalar fields 
only gravitationally. In a recent paper Coley and van den Hoogen qualitatively analyzed the autonomous system 
of two scalar fields of the kind discussed above in a curved FLRW universe and showed that the system has an 
equilibrium point compatible with a stable phase of power-law inflation. This feature persists even if a perfect fluid 
with baryotropic equation of state is allowed in the picture 0| . Again the authors assumed that the fields do not 
interact among themselves nor with the matter fluid. 

At first sight the presence of a matter fluid (such as a sea of relativistic particles) may seen as rather irrelevant 
because the fast inflationary expansion will very soon dilute away these particles. However, as shown by Berera |^ 
if a coupling between the inflaton fleld and the matter fluid is assumed, things change drastically to the point that 
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when inflation ends the Universe is far from being so cool as predicted by ordinary inflation models. In fact it can be 
hot enough to resume the radiation dominated phase of standard big-bang cosmology, thus the reheating phase is no 
longer necessary. 

The main target of this paper is to study the scenario of warm inflation with n scalar fields having exponential 
potentials and interacting with the matter fluid. Section II considers first the more simple case of a single scalar 
field plus a dissipative matter fluid, Then the analysis is extended to n scalar fields. In section III the scalar fields 
are allowed to decay into the matter fluid thereby the temperature of the latter does not fall drastically and so the 
warm inflationary phase can be followed by the conventional radiation dominated period without any intermediate 
reheating. Section IV studies the perturbations in energy density and entropy brought about by warm inflation; they 
do not conflict with the observed temperature anisotropics of the cosmic microwave background radiation. Finally, 
section V summarizes our main findings. Units have been chosen so that c = SttG = 1. 



II. THE SYNERGISTIC MECHANISM 

Let us assume a FLRW universe filled with a self-interacting scalar field plus a dissipative matter fluid. The 
stress-energy tensor of this mixture is 

Tab= {p+P + TT)UaUiy+ [p + TT)gab, (1) 

where p — p„i + and p — Pm + P<f>- Here pm and p„j are the energy density and pressure of the matter fluid with 
equation of state given by Pm — {jm — ^)Pm and with its baryotropic index in the interval 1 < 7^ < 2. Likewise 
and p^, the energy density and pressure of the minimally coupled self-interacting field (f), i.e., 

P^ - l<i^' + V{(b) , p^ = ~ , (2) 

are related by an equation of state similar to that of the matter, namely p^ = (7^ — so that its baryotropic 

index is given by 

where for non-negative potentials V{(f>) one has < 7^ < 2, and an overdot means derivative with respect to cosmic 
time. In general 7^ varies as the Universe expands, and the same is true for 7,„ since the massive and massless 
components of the matter fluid redshift at different rates. 

The Friedmann equation together with the energy conservation of the normal matter fluid with bulk dissipative 
pressure and Klein-Gordon equation can be written as 

n,„. + n^ + nK = 1, (i^ = i,o,-i), (4) 

Pm + ( + — ] P,n = 0, (5) 

V PrnJ 

p^ + SH"f^p^ = 0, (6) 

where H = a/a denotes the Hubble factor. We have further introduced the density parameters ^Im = Pm/Pc 
f20,= P(t,/pc, with pc = 3i?2 the critical density, and = —K/[aH)'^. 

In terms of these quantities we can introduce an overall baryotropic index 7 

-fVi = 7,„ri,„ + 70rj0 , (7) 

where we have made use of the definition f2 = + 17^. The flatness problem is solved by the attractor solution f2 = 1 
of equation (|^). In addition, the ratio 17^/17™ becomes asymptotically a constant, meaning that the matter content 
of the universe does not dilutes altogether as inflation proceeds. In it was shown that these constant solutions are 
stable in the (O, Q.m, ^<f,) space. The fixed point solution Vl — 1, Q.m = ^mo and Jl^ = fi^o, respectively, of equations 
(@)"® is obtained when the partial adiabatic indices and the dissipative pressure are related by 
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7m + — =70 = -X77^, (8) 

accordingly the smaller 7^, the larger the dissipative effects. 

Typically bulk viscosity arises in mixtures either of different particles species, as in a radiative fluid, or of the 
same species but with different energies, as in a Maxwell-Boltzmann gas. Physically, we can think of tt as the 
internal "friction" that sets in due to the different cooling rates in the expanding mixture. Any dissipation in exact 
FLRW universes has to be scalar in nature, and in principle it may be modelled as a bulk viscosity effect within a 
nonequilibrium thermodynamic theory such as the Israel-Stewart's pO| |. In that formulation and under certain general 
circumstances, the evolution equation for the bulk dissipative pressure takes the form 

TT + TTT = -i^H , (9) 

where the positive-definite quantity ^ stands for the phenomenological coefficient of bulk viscosity, T the temperature 
of the fluid, and r the relaxation time associated to the dissipative pressure -i.e., the time the system would take 
to reach the thermodynamic equilibrium state if the velocity divergence were suddenly turned off |l^]. Usually ^ is 
given by the kinetic theory of gases or a fluctuation-dissipation theorem or both ||l^ . Expression (H) has been widely 
used in the literature p^ ] and it meets the requirements of causality and stability to be demanded to any physically 
acceptable transport equation fl^ . 

Combining (||) and we obtain the equation of motion of the attractor solutions satisfying flatness, accelerated 
expansion and the non-dilution condition 

.-(| + 3.„h)+H + 2|.H»-^H = 0. (10) 

Here v = (tH) ^ is the number of relaxation times in a Hubble time - for a quasistatic expansion ly is proportional 
to the number of particle interactions in a Hubble time. Perfect fluid behavior occurs in the limit ;/ 00, and a 
consistent hydrodynamical description of the fluids requires v > 1. 

The problem of a homogeneous scalar field driven by a exponential potential 

V^ = ^oe-^^ (11) 
minimally coupled to gravity in a flat FLRW spacetime with a linear viscosity coefficient 

e-Coi?, (12) 

where Vq and ^0 are constants, has the solution 

a = ao{t/toT (13) 
= 0oln(tAo). (14) 

The quantities 

(15) 



(1 ^ma) 

Co = ^ma (7m " l4>a) [l " 37^a , (16) 



and 



70a = ^ (17) 

are obtained by solving the system (^-(||), ( p^ and ( pi] ) -the subindex a stands for asymptotic value of the 
correspondin g qu antity. The power-law expansion (|13|) will be inflationary for < 2(1 — il„ia)^^. 
Rewriting ( |lO|) in terms of the field baryotropic index 7^, we get 



70 = 370 -{1^ + 37m)70 + (im - ^ 



(18) 
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where a prime indicates derivative with respect to r] = In a. When the phenomenological coefficient of bulk viscosity is 
given by ( |l2|) and ([l6|), that is ^a, Eq- ( |l8| ) admits the constant solution 7^ = 70a- It gives an accelerated expansion 
in the late time regime when 7^^ <2/3. As^>0 and "fm ^ 1, the hydrodynamical parameter v is restricted to 
u > 370Q. The case of constant 1^0 arises for instance in a radiating fluid, and the nearly linear regime, with slowly 
varying v and 7„j, was already investigated in the quasiperfect limit, corresponding to — > 
To analyze the stability of the solution 7^ — we insert (^2|) in (|l^) to obtain 

7; = 3 (7^ - 70a) -{'^ + 37m) (70 - 70a) (19) 

As 70a < 2/3, V > max(370a, 1) and 7,„ > 1, Eq. ( p^ shows that d^^/djif, < in a neighborhood of 70a. Hence this 
constant solution is asymptotically stable, showing that all solutions of Eq. (|lo|), that is all the accelerated attractors 
of the system (^), (||), (|^, (||), are themselves attracted towards the constant solution 70 = 70a provided that they 
satisfy ~ Ca at late time. 

Let us now assume that instead of having just one scalar field, we have n homogeneous non-interacting scalar fields, 
with exponential potentials Vi = Voi e"^'*^'. In that case the Einstein-Klein-Gordon equations can be recast as 

n 

n^ + J2n^, + nK^i, = 1,0,-1), (20) 



Pm + 3H (j„i + —) P^ = 0, (21) 

P0. + 3i770^ p0i = (i = l,2,...n). (22) 

As it stands the problem in its full generality is rather involved. Therefore to bring it to a form amenable to 
analytical treatment we shall consider henceforth the a simplified version characterized by Ai = ^2 = • • • = = A 
and Vbi = V02 = • • • = Vbn = Vb . As we shall see in the next section, we can expect in this particular case that all scalar 
fields share the same asymptotic limit. So, in the remaining of this section we will assume (pi ^ (j)2 ^ ■ ■ ■ — (pn = 
so that Vi = V2 — ■ ■ ■ = Vn = V = Vq e""*"^ and equation ( pO| ) becomes into 

n„, + nil^^ + = I, {K = 1,0,-1), (23) 



while (E2h simplifies to m). Following parallel steps to that leading to Eqs. 



it can be seen that the Einstein- 



Klein-Gordon system has the spatially flat power-law attractor solution ([13) but now with 

2n 



(1 ^ma) ^ 



(24) 



showing that the n scalar fields p cooperate to a stronger inflation. 

Let us now assume that the n homogeneous scalar fields, (pi are driven by a general potential V — V{(j)i). In that 
case the Einstein-Klein-Gordon equations are 

3^^ = ^E0? + ^ + Pm-^. (25) 

4>, + m^, + 1/0, = 0. (26) 

along with equation ( ^l|) (1^0^ stand for dV/dpi). From these equations we get 

i^=-^f]'?^»'-^(7mP™ + 7r) + J. (27) 
In order to investigate the stable scalar field configurations it is expedient to introduce the ancillary quantity 

n9a 
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which reduces to w = 1 for the completely symmetric configuration cjji 
differential equation for uj: 



Using (E§-(28) we find the 



(jj 



(29) 



(In this section no summation convention applies to repeated Greek indices) . If we further assume that the potential 
satisfies the condition 



V ^nV 



0„ 9a, 



then equation (|29|) becomes 



^^2^(w-l). 



(30) 



(31) 



This has a fixed point solution, namely oj = 1. Further; with the aid of (^)-(^) the general solution of (|3l| ) can be 
found in terms of the scale factor, the matter fluid and the bulk viscous pressure 



nc 



(jj 



2a6 



(32) 



where c is an arbitrary integration constant. Evaluating ( |32[ ) in the asymptotic regime of the attractor solution 
70 — lij}a and the potentials Vi — Voie~^^^% which asymptotically satisfies the condition (^0|), it can be easily shown 
that the particular solution w = 1 is an attractor for evolutions that behave asymptotically as a oc with a > 1/3. 
On the other hand this result strongly suggests that the special case in which all scalar fields are equal may be the 
late-time attractor of more general scenarios. 



III. WARM INFLATION 



Warm inflation arises when a strong enough coupling between the scalar field and matter fluid (which we shall 
assume perfect) exists. The former decays into the latter (which acts as a thermal bath) while the inflaton slowly 
rolls down the potential . The decay is phenomenologically implemented by inserting a (usually constant) friction 
term F in the equation of evolution for 

+ 3(ff + F)0 + F'(0) = 0. (33) 

We adopt this picture except that (i) we consider no slow-roll (although it can be straightforwardly incorporated), 
and (ii) rather than a single field we have n scalar fields all of them with identical exponential potential. Here 
stands for any of these fields. Accordingly ( p3| ) is not just a single equation but n identical equations; besides for 
mathematical simplicity we assume the same friction term for each field. 



Obviously the coupling between the n scalar fields and the matter fluid introduces a source term in the energy balance 
equation for the latter 

Pm + i^rnHpm = 3F0^. (34) 

So there is a continuous transfer of energy from the scalar flelds to the matter adjusted in such a way that the 
former experience a damped evolution and give rise to a nearly isothermal expansion. Accordingly, like in standard 
warm inflation, no reheating mechanism is needed at the end of inflation. Moreover, thermal rather than quantum 



fluctuations produce the primordial spectrum of density perturbations p9-23|. 

We can identify the phenomenological coupling with an effective dissipative pressure tt* along the evolution on the 
attractor. Then comparing (^ with ( |3^ ) we get 

TT = —. (35) 
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In the case at hand the attractor condition (g) becomes dynamic because the starred magnitudes include the interaction 
between the scalar field and radiation 



(36) 



where 7^ = 7^ — {-k* / p^). Using this together with (|^) and (H) it follows that 



(37) 



(bear in mind that 7m > 7^. This implies that the scalar field evolves with an "effective" expansion rate H = {1+R)H. 
Once converged to the attractor solution, R becomes a constant and the effective power-law exponent results larger 
by a factor (1 + i?): 



2n (1 + R) 

(1 ^ma) ^ 



(38) 



This means that this interaction between the scalar fields and radiation further assist to infiation. 

Let us investigate the entropy production S along this era. Using (|3^) and (||) we get that matter redshifts 
as pra = Pmeifll a-e)~'^'^'^ (subiudcx 6 mcans evaluation at inflation exit). As it is customary in warm inflationary 
scenarios we assume that this matter behaves as radiation pm = {tt^ /30)gT'^, where g is the effective number of 
relativistic degrees of freedom, so we have 



Then, using (§, (|3|), (|l|) and S = n*^/^T, we get 

S{t) = S, 



(39) 



(40) 



where 



Tl 



16 



-4/3)' 



1/2 



2/3 



(41) 



is the characteristic time for the entropy density to attain a constant value Se and is directly related to the time for 
the start of warm inflation. 

We see from ( p9| ) that the requirement of a nearly isothermal expansion along the inflationary era imposes a 
constraint on the value of 7^. Namely, if Ti is the temperature of the radiation fluid at the beginning of inflation, Tg 
at the end, and N is the number of e-folds, we have 



(42) 



Assuming TijT^. ~ 0(1) we find that 7J ~ 10^^. From ( |l7|) and (|24| ) we see that this can be easily achieved, with n 
in the range of 10 — 100, without fine-tuning the potential slope. 

Inflationary scenarios need to achieve a graceful exit from inflation. In our case this is not a problem since the ratio 
P4>/ Pm is not exactly constant, as can been seen from ( |33| ) and (p4[). Therefore, the continuous transfer of energy from 
the n scalar fields to the matter fluid slowly increases the energy of the latter and decreases that of the former. Thus 
the acceleration equation a/a= —^[p + 3p) implies that the universe ceases to accelerate when both energy densities 
equalize. This criterion for ending inflation coincides with Taylor and Berera's [M. 
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IV. EVOLUTION OF PERTURBATIONS 



This section considers the evolution of energy density, entropy and curvature fluctuations in the perturbative long- 
wavelength regime during the attractor era. Scalar perturbations are covariantly and gauge-invariantly characterized 
by the spatial gradients of scalars. Energy density inhomogeneities are described by the comoving fractional density 
gradient 

= ^ , (43) 
P 

where stands for the covariant spatial derivative T)jAi... — hj'^h/ ■ ■ ■ VfeA;.... The scalar part S = aD^Si = {aD)"^ p/ p 
encodes the total scalar contribution to energy density inhomogeneities. It relates to the usual gauge-invariant density 
perturbation scalar e„i through 5 = V^Em, where is the Laplacian for the metric of the 3-surfaces of constant 
curvature [^5y26[ . Also the comoving expansion gradient, the normalized pressure gradient, and normalized entropy 
gradient are defined by [ p4|p^ 

^ ^ aDiP anTDiS ,,,, 
e,^aD,e, p, = ^, e, = (44) 

P P 

n being the particle number density, T the temperature, and s the specific entropy per particle. The evolution equation 
for scalar density perturbations reads 

J + (8 - 67 + Zcl) 5 - j-lO + 147 - 87^ - &cl+ 



(1 - 3 (7 - 1)' + 2c2] fc} 5 - c^D^J = S[e] + S[^*] + S[q] + S[a] , (45) 



where 

dp ) ^ nT \ds 

are, respectively, the adiabatic speed of sound and a non-baryotropic index. The sources in the right-hand side of Eq. 
( p5| ) arising, respectively, from entropy perturbations, bulk viscous stress, energy flux, and shear viscous stress are 
given in p7| . Since in our case there are no shear viscous stress (ct^ ~ 0) and S[q\ vanishes by choosing the energy 
frame (qi = 0), we reproduce here only the expressions for S[e] and S[7r*]: 

S[e] = r (Si^i/^ +d2) e, (47) 
S[n*]^-{?,KH^ + \i'^)B, (48) 



where the scalar entropy perturbation 



and the dimensionless perturbation scalar 



e = aD'e^ = D^s (49) 



(50) 

P 

related to the inhomogeneous part of the bulk viscous stress, were defined. Also, the entropy perturbation equation 
in the energy frame is 

e + 3ff (c^ -7 + 1 + r) e = -3i/S . (51) 

The coupled system that governs scalar dissipative perturbations in the general case is given by the energy density 
perturbation equation (^), the entropy perturbation equation (|5l|), the equation for the scalar bulk viscosity (|9|), and 
the equation for temperature perturbations. 

When only bulk viscous stress dissipation is present, the coupled system can be reduced to a pair of coupled 
equations in S (third order in time) and e (second order in time). For a flat background, the equations are p7|]: 
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T 'c^' + [1 + 3 (2 - 7 + c^) rff] ^ + {s - 67 + 3c2 + 3r (c^)' 



-i [-14 + 757 - 487^ + (2I7 - 30)c2] tH} S 



H'^ |6 - IO7 + 57^ - - 4r (c^) ' 



-2 [-6 + I87 - 157^ + 57^ + (6 - 287 + 107^) c^] tH} S 



P7 V y p7 



(l-37ri/)c.^ + r(c.^)-+3(|) 



and 



(1 - 37Ti7) r + Tf + 3 ( 1^ 



D^e. 



re + [1 - I (-2 + 37 - 2c^ - 2r) tH] e 



(52) 



-3i? 



7 — 1 — Cg — r + 37 (7 — Cg) ri/ + 7" (cg + r) 



7 \9s 



e = d H 

7 7 



(53) 



We shall consider here the evolution of the energy density and entropy perturbations in the attractor stage with 
the conditions r — and dv/ds — 0. Together with Eq. (12) they imply 



dp) s \dPrn 



(54) 
(55) 



In this case Eq. ( |5^ ) decouples to give 

Cl v: , C2 1 C3 



t t-^ 



W5 



and Eq. (53) becomes 



.. Cs . Cg Clo j Cll c 



(56) 



(57) 



where the constant coefficients ci . . . cn depend upon the parameters of the model: v, ki, K2, Vim, a, 7, Cg , and the 
value of the scale factor ae at the exit from inflation. For our purposes, only the explicit expression for ci, cs, and cg 
are relevant, being 
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Ci — av + 2ia (2 — 7 — Cg) , 



j.--(37-2-2c2) 



Cg 



-3«2 [j. (7 - 1 - c^) + 37 (7 - cl)] 



(58) 



We deal with the system (|56|),(57) by performing separation of variables in the form 6 = 6x 5t and e = ex Cf , where 
5x and ex depend upon the spatial variables while 5t and et are functions of the coordinate time t. Then, Eq. ( [5^ ) 
can be recasted as 



5t 



2^ 



C_6 



— + C7 



(59) 



which can only hold if 



(D^ - /i) (5, = 



(60) 



PL being an arbitrary constant. In Fourier transformed space it becomes an identity that holds for an arbitrary 
amplitude 5k with /i — —k"^ /o?, where k/a is the physical wavenumber. Also, Eq. (^ ) leads to 



(et/et) + cs t {et/et) + cg 



cio 



[St/e^ + cii ((5t/et) 



(61) 



which requires e^, = A (5^; , with A a constant. Then the evolution equation for mode /i becomes 



Cl 



t 



C2 



i2 



C5-M(f + C7) 



(5t = 



(62) 



As we are interested in the asymptotic behavior of the perturbations along the attractor regime, it suffices to 
consider the dominant terms in ( |62[ ) for large time. As a > we have 



"5t +j5t - Ai' t'" C4 5t - c5St = 



(63) 



To study the asymptotic evolution of the long-wavelength modes much larger than the Hubble scale {k/aH <^ 1), we 

expand 8t in powers of /i^ as 8t = <5j°^ + jj? Then, replacing this expression in Eq. (|63| ) and retaining terms up 
to first order in /i^, we obtain 



t 



(0) 



Its zeroth-order solution is 



(l-ci)(2-ci) 



,2-ci 



A2t + A3, 



(64) 



(65) 



for Cl ^ {1,2} and being Ai , i = 1,2,3 arbitrary integration constants. Then, S^^^ satisfies the inhomogeneous 
equation 



(66) 
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whose general solution has the form 



c(l) 



,2-ci 



Bo t + B, 



(l-ci)(2-ci) 
+a2 + aa t^^+^-^i + 



ait 



2a+4-ci 



(67) 



where the coefBcients i = 1, 2, 3 are integration constants and a^, i = 1, . . . , 4 depend on a, ci, and the Bi. Bearing 
in mind that a ^ 1, we find that the dominant mode of the energy density perturbations in the super-Hubble regime 



grows asymptotically like i2a+4 



oc a 



2+(4/a) 



a , independently of the wavenumber, until inflation exit. 



On the other hand, Eq. (53) becomes, in the leading regime, 



Cg . Cg 

et + -e, + -e, 



ci2r 



1+- 7-1 



2il 



whose solution is 



et{t) 



B.r 



Cut 



2a+3 



7-1 



(68) 



(69) 



where A1.2 are the roots of the equation + (cg — 1)A + C9 — 0, -B4, B^ are arbitrary integration constants, and C12, C13 
are functions of the parameters and the previously defined integration constants. It follows that the dominant mode 



of the entropy perturbations also grows as a 



2+(4/a) 



for 1^ = 0(1). 



To deal with the evolution of the curvature perturbations generated by these energy density and entropy fluctuations 
we will turn to the standard metric-based gauge-invariant approach. In a particular choice for slicing of space-time 
named the longitudinal gauge, the metric describing the inhomogeneous perturbations of the spatially flat FLRW 
background takes the simple form 



ds^ = (1 + 2$) dt^ - a^(t)(l - 2*)(5.y dx' dx' , 



(70) 



in terms of the gauge invariant Bardeen potentials and ||2^;|2^. Besides, when the shear stress vanishes, it follows 
from the equations of motion for the gauge invariant variables that $ = . So, just a single scalar degree of freedom, 
say is required to describe linear perturbations of the metric. We get two second order equations, namely 



2 ^' 



(71) 



+ 3aH (1 + cl) $' 



[2a {aH)' + (1 + 3c^,)a'^ H"^] $ 



-a^ {Sp + Stt* - cl6p) , 



(72) 



where the prime denotes derivative with respect to conformal time, and the source terms contain the energy density 
perturbation 5p, the equilibrium pressure perturbation dp and the dissipative pressure perturbation Jtt*. 

When the source term of Eq. ( fz^ ) vanishes, and the scales are larger than the Hubble radius such that the 
spatial gradients can be neglected, Eq. ( [7^ ) can be recast in terms of the curvature perturbation on uniform-density 
hypersurfaces p^~|3l|| 



C = $ + i? 



5p 



2i/-i$ + $ 
3 7 + 7r*/p 



(73) 



as a conservation law C = ||3^,Q. We consider that the equilibrium pressure perturbation is isentropic during 
warm inflationary and radiation dominated eras, and that the dissipative pressure perturbation switches off during 
the transition between both eras with a relaxation time that is a fraction of a Hubble time. Hence, one should expect 
that C, soon vanishes along this transition so that the value Ce at inflation exit on super-Hubble scales may be equated 
to that at reentry of long wavelength modes to Hubble radius during the radiation- or matter-dominated eras. To 
find this value we need the evolution of <& along the attractor era for long-wavelength modes. Again, neglecting the 



spatial derivative term and using the relationship e,„fc = —StSk/k 
becomes 



Sk/k , the Fourier transform of Eq. (|7l|) 
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Its solution 



shows that <& grows asymptotically as o? during the attractor regime. Inserting this result back in ( |7^ ) we find 
that C a-lso grows as a? during this regime. This nonconservation of the super-Hubble curvature perturbations is a 
consequence of the growth of the entropy perturbations, which, in its turn, is due to dissipation effects as Eq.([5l|) 
shows. Then, taking into account that ^ 1, we find that C, — 3a$ ~ —aem/'^ at inflation exit. Hence the 
power spectra of C and $ at that moment are proportional to the power spectrum of the primordial energy density 
perturbations. 

For a mode that crosses outside the Hubble radius at scale factor aA during inflation and reenters to the Hubble 
radius at scale factor as during the radiation dominated era, we find the number of e-foldings before the end of 
infiation Na = a/{a. — 1) ln(aB/ae) by using the continuity of the energy density at the turnover between the warm 
inflation and the radiation dominated eras. Then, for this mode of perturbations, the regime of growth as starts at 
scale factor oi 3> oa, when both aH 3> k and the evolution is close to the attractor, and it continues until inflation 
exit. Then we obtain 



q/(q-1) 



Ob 



ai 



a \(3aj\j J l3aaN 



(76) 



where a — ai/a^, a^r is the scale factor at the start of nucleosynthesis era, corresponding to a temperature Tjv — 
IMeV, (3 — a^/aM and we are taking a 3> 1. As in this warm inflationary scenario there is no need to accommodate a 
reheating stage, inflation may end shortly before nucleosynthesis and we may take safely (3 ~ 10~^. Then we obtain an 
upper bound on the growth of perturbation modes crossing inside the Hubble radius during the radiation dominated 
era by considering the matter-radiation equality scale ~ lOOMpc, corresponding to a cluster of galaxies. As 
o-eq/o-N — 10^, taking cr >,10'^ we find that eme/^mi = Ce/Ci^lO"'^^- Similarly, for a mode crossing inside the Hubble 
radius during the matter dominated era, we obtain 



ai ~ 4(3(7 \aN ) \ aeq 



1/2 

(77) 



Thus, scales fc^J^ ~ lO'^Mpc, corresponding to the observable universe, give the upper bound on the growth of 
perturbation modes Ce/Ci~10^'^- 

After the inflationary stage, we recover the standard picture of conserved isentropic curvature perturbations and 
well known calculations show that the density contrast at Hubble scale entry {Sp/ p)k=aH is proportional to C, or 
equivalently to the comoving curvature perturbation TZ, with a proportionality factor of 2/5 during presureless matter 
dominated era, or 4/9 during the radiation dominated era ]3^ . So, using COBE normalization for the power spectrum 

of curvature {2/5)T>^'^ = 1.91 x 10"^ at the scale fc-^ ~ lO^Mpc, we find that at inflation exit emk = 9.55 x 10 ^/a ~ 
lO^*' for a ~ 100. Besides, due to the proportionality of power spectra, the observed bound on the spectral slope of 
curvature perturbations n = 1.0±0.2 implies the same bound on the spectral slope of the energy density perturbations 
at inflation exit. Recalling that the amplitude Sk is an arbitrary function, these observational bounds impose no fine 
tuning constraints whatsoever on the parameters of the scalar field potential. 



V. DISCUSSION 



We have proposed a new inflationary scenario whose main ingredients are n scalar fields and a dissipative matter 
fluid. The former decay into the latter at a high rate F. While no single scalar field can achieve inflation by its own 
they cooperate synergistically to produce it. We have derived the attractor condition Eq. (H) and shown that the 
presences of dissipation does not spoil the linear relationship between the power-law exponent a and the number of 
fields, preserving the stability of the symmetric configuration of n identical fields. 

We have described the interaction between the scalar fields and the radiation fluid in the warm inflationary scenario 
by means of an effective bulk dissipative pressure and generalized the attractor condition. Likewise, we have resorted 
to the synergistic mechanism to calculate the production of entropy and the evolution of temperature. The exit 
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temperature results lower but approximately of the same order than the initial temperature. This renders the reheating 
phase redundant. 

Further, we have found that the combination of the synergistic mechanism and the decay of the scalar fields 
into the matter fluid produces significant entropy perturbations with proportional spectral amplitude and dominant 
mode evolution to that of energy density perturbations on large wavelength scales until inflation exit. This steep 
growth contrasts with other models of inflation where long wavelength curvature modes evolve isentropically (see e.g. 
psf , Jsst); however in our case there is a continuous transfer of energy from the scalar fields to matter. Moreover, 
observational bounds on the curvature perturbations at Hubble scale entry do not force on our model any slow-roll 
constraints on the scalar field. 
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